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ABSTRACT 

We present the results of a programme of thermal-infrared imaging of nineteen 2: ~ 1 radio 
galaxies from the SCR and 3CRR (LRL) samples. We detect emission at L' (3.8 /im) from 
four objects; in each case the emission is unresolved at 1-arcsec resolution. Fifteen radio 
galaxies remain undetected to sensitive limits of L' ~ 15.5. Using these data in tandem 
with archived HST data and near-infrared spectroscopy we show that three of the detected 
'radio galaxies' (3C 22, 3C 41, and 3C 65) harbour quasars reddened by Ay ^ 5. Correcting 
for this reddening 3C 22 and 3C 41 are very similar to coeval 3C quasars, whilst 3C 65 
seems unusually underluminous. The fourth radio galaxy detection (3C 265) is a more highly 
obscured {Av ~ 15) but otherwise typical quasar which previously has been evident only in 
scattered light. We determine the fraction of dust-reddened quasars at 2; ~ 1 to be 28^13% 
at 90% confidence. On the assumption that the undetected radio galaxies harbour quasars 
similar to those in 3C 22, 3C 41 and 3C 265 (as seems reasonable given their similar narrow 
emission line luminosities) we deduce extinctions of Ay ^ 15 towards their nuclei. The 
contributions of reddened quasar nuclei to the total K-hand light ranges from ~ per cent 
for the non-detections, through ~ 10 per cent for 3C 265 to ~ 80 per cent for 3C 22 and 3C 41. 
Correcting for these effects does not remove the previously reported differences between the 
K magnitudes of 3C and 6C radio galaxies, so contamination by reddened quasar nuclei is 
not a serious problem for drawing cosmological conclusions from the K-z relation for radio 
galaxies. We discuss these results in the context of the 'receding torus' model which predicts 
a small fraction of lightly-reddened quasars in samples of high radio luminosity sources. We 
also examine the likely future importance of thermal-infrared imaging in the study of distant 
powerful radio sources. 
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1 INTRODUCTION 



Luminous extragalactic radio sources are associated with ei- 
ther quasars or radio galaxies, the latter class of object lack- 
ing the bright non-stellar continuum and prominent broad 
emission lines that characterize the former. For many years 
astronomers have speculated that at least some radio galax- 
ies harbour an active nucleus which, although obscured from 
view along the line of sight, would be observed as a quasar if 
seen from certain other, favourable directions (e.g., Scheuer 
1987). There is now overwhelming observational evidence 
that this is the case. Most notably, spectropolarimetry has 
revealed the presence of quasar-like broad emission lines in 
the polarized flux spectra of a number of radio galaxies (e.g., 
Tran et al. 1998 and references therein). In such cases the 
broad lines are hidden from direct view but are scattered 
into our line of sight by a screen of either dust or electrons. 
However, the more general proposition that all radio galax- 



ies harbour obscured quasar nuclei has not yet been con- 
firmed. 

It is generally believed that the unobscured lines of sight 
lie close to the axis of the twin radio jets, and that quasars 
are therefore seen more nearly 'pole-on' than their radio 
galaxy counterparts (e.g., Barthel 1989). This is supported 
by the presence of bright radio cores and one-sided jets in 
quasars that are weak or absent in radio galaxies (e.g., Owen 
& Puschell 1984), a result which is most naturally explained 
by Doppler boosting of emission from relativistic jets mov- 
ing close to the line of sight. Even more compelling evidence 
comes from optical polarization studies that have revealed 
many radio galaxies in which the rest-frame ultraviolet light 
is strongly polarized with its electric vector lying perpendic- 
ular to the radio axis (e.g., Tadhunter et al. 1992). These 
observations suggest that quasar Hght can emerge relatively 
unhindered along the radio axis, and reach the observer af- 
ter scattering, whilst direct transmitted radiation from the 
quasar nucleus is effectively blocked. Such processes are also 
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at least partly responsible for the so-called 'alignment ef- 
fect' (Chambers, Miley & van Breugel 1987; McCarthy et 
al. 1987) — the tendency for the elongation of optical and 
radio structures along a common axis. 

In an orientation-based unification scheme such as this, 
the obscuring material lies preferentially in the plane per- 
pendicular to the radio jets, and a central parameter is the 
angle between the radio axis and the line of sight which 
grazes the edge of the obscuring material (generally referred 
to as the 'torus'). A simple model which is often considered 
relates this critical angle 9^ to the inner radius of the torus, 
r, and its half-height, h, via tanSc = r/h (see, e.g., Simp- 
son 1998). In this model there are physical reasons why r, 
and hence dc, might be expected to correlate with quasar 
luminosity: if the value of r is determined by the radius at 
which dust sublimates in the radiation field of the quasar 
then it should increase with the photoionizing luminosity of 
the quasar as Lpifot (Lawrence 1991). This picture where the 
torus opens up as the luminosity of the quasar increases in 
known as the 'receding torus' model. Strong links between 
the intrinsic UV/optical luminosity of the quasar nuclei of 
radio sources and extended radio luminosity are now firmly 
established (Rawlings & Saunders 1991; Falcke, Malkan & 
Biermann 1995; Serjeant et al. 1998; Willott et al. 1998c) 
in the sense that Lphot ex: Z/17I where Lirg is the 178-MHz 
(extended) radio luminosity. It therefore seems highly plau- 
sible that r and hence 6c will rise, albeit slowly, with radio 
luminosity (tan 6*0 oc iiTs), 

The brightest (3C) radio sources at z ~ 1 are about 100 
times more radio luminous than 3C sources at z '^ 0.1, so 
assuming h to be independent of Lits and z, we expect r and 
hence tan 9c to be a factor ~ 4 greater for these objects. The 
quasar fraction should therefore be correspondingly larger, 
and, since lines of sight grazing the torus will subtend a much 
lower solid angle, the fraction of lightly-reddened quasars 
should be much reduced (see Hill et al. 1996). Based on their 
detection of a high fraction of lightly-reddened quasar nu- 
clei from near-infrared spectroscopy of Paa in low-redshift 
(0.1 < z < 0.2) 3C radio sources, HiU, Goodrich & De- 
Poy (1996) found 9c ~ 20°, and used the receding torus 
model to make quantitative predictions about the distribu- 
tion of nuclear extinctions in higher-redshift, more radio lu- 
minous, sources. Their predictions were in rough agreement 
with Economou et al.'s (1995) discovery of only one broad 
Ha line in a preliminary analysis of a sample of ten radio 
galaxies at z '^ 1. However, as Rawlings et al. (1995) dis- 
cuss, spectroscopic methods of detecting hidden quasar nu- 
clei lack the sensitivity to provide a good general method of 
characterizing the angular distribution of reddening towards 
quasar nuclei in radio galaxies at 2 ~ 1. This is because of 
the large optical depth expected towards rest-frame optical 
lines like Ha, and the poor sensitivity to weak rest-frame 
near-infrared lines like Paa, particularly once they become 
redshifted beyond the iC'-band (at z > 0.3 in the case of 
Paa). A definitive test of the receding torus model for high- 
redshift objects has therefore yet to be performed. 

The distribution of reddenings to the quasar nuclei may 
also have important cosmological implications, since the 
near-infrared Hubble diagram for radio galaxies (Fig. y) — 
and the common assumption that the near-infrared emission 
of radio galaxies is dominated by starlight — has been used 
to infer constraints on the formation and evolution of ellipti- 



O 
C\2 



^ 



LO 




Figure 1. The near-infrared Hubble Diagram for radio sources. 
The objects are taken from the 3CRR (LRL) sample (squares; 
Laing, Riley &; Longair 1983), and the 6C sample (triangles; 
Eales 1985, slightly revised by Rawlings et al. 1998b); quasars 
are represented by filled symbols, and other objects by open sym- 
bols. We define quasars are having a nuclear point source with 
Mb < —23; with this criterion, all quasars display broad optical 
emission lines, but some of the other objects, which we term ra- 
dio galaxies, also display broad lines, either weakly or in polarized 
fiux only. X-band photometry for the LRL sample has been taken 
from Best et al. (1998), Lilly & Longair (1982, 1984), Lilly, Lon- 
gair & Miller (1985) and Rawlings et al. (1995). Although 73 of 
the 94 2 > 0.1 LRL radio galaxies have K-ha.nd photometry (and 
in all but one case the lack of _fsr-band data stems from the in- 
ability to observe high-declination objects at UKIRT), almost all 
the quasars lack such data; to plot the quasars we have estimated 
K-haad magnitudes using optical photometry from LRL, and an 
extrapolation based on the canonical quasar spectrum described 
in Section hi Note that contributions from the host galaxies of 
quasars, and the possibility of their suff'ering significant reddening 
mean that these estimated K values should be treated as upper 
limits, i^-band photometry and redshifts for 6C sources are from 
Eales et al. (1997) and Rawlings et al. (1998b) respectively. 39 
of the 41 z > 0.1 6C radio galaxies in the 6C sample have K- 
band photometry; we again estimate if -band magnitudes for the 
7 quasars using optical photometry (Eales 1985) and the canon- 
ical quasar spectrum. The lines show the loci of a model quasar 
with Mg = —26 reddened by Ay = 0, 2, 7, 15; the lowest curve is 
for Av = 0, increasingly higher curves represent increasing values 
of Ay. A few objects of special interest are labelled. 
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cal galaxies (e.g., Lilly 1989; Eales et al. 1993). The discovery 
that the K-hanA magnitudes of 6C radio galaxies, i.e., those 
selected at a ~ 6 times lower radio flux limit than 3C, are 
~ 0.6 magnitudes fainter than those of 3C radio galaxies at 
the same redshift (Eales et al. 1997) suggests that there is 
some close link between radio luminosity and near-infrared 
luminosity. Either, stellar luminosity — and by inference 
stellar mass — correlates with radio luminosity, or there is 
some important component of the near-infrared luminosity 
which can be ascribed to a hidden quasar nucleus. Recent 
discussions of these possibilities can be found in Eales et al. 
(1997), Best, Longair & Rottgering (1998) and Roche, Eales 
& Rawlings (1998). Although there is direct observational 
evidence from near-infrared polarimetry (Leyshon & Eales 
1998) that scattered light makes only a small (few per cent) 
contribution to the near-infrared luminosity of radio galax- 
ies even when it makes a very significant contribution in 
the ultraviolet and optical, it has not yet been conclusively 
shown that reddened quasar light is also a negligible compo- 
nent. As illustrated by Fig. hi, such quasar nuclei would need 
to be subject to only a fairly limited amount of reddening 
{Av ~ 2-7) in order not to possess obvious quasar features 
in their optical spectra whilst still contributing to the K- 
band luminosity at the level of radio galaxies lying close to 
the mean K-z relation for radio galaxies. The existence of a 
substantial population of moderately-reddened quasars has 
been suggested by Webster et al. (1995), although both the 
reality and importance of this population have been ques- 
tioned (Serjeant & Rawlings 1996; Benn et al. 1998). 

In this paper we report on the results from a project 
designed with two principal aims. The first aim was to de- 
termine the contribution of transmitted quasar light to the 
K-z relation for z ~ 1 radio galaxies, and hence determined 
whether reddened quasar nuclei can explain the difference in 
near-infrared luminosities between 3C and 6C radio galax- 
ies. The second aim was to constrain the range of Ay to- 
wards the quasar nuclei of luminous radio sources at 2 ~ 1, 
and to compare this range with the predictions of the re- 
ceding torus model. The project was based around thermal- 
infrared imaging since, as discussed by Rawlings et al. (1995) 
this is for these purposes a much more powerful technique 
than spectroscopy. Its utility rests on the dramatic diver- 
gence of the rest-frame spectral energy distributions (SEDs) 
of starlight and quasar light redward of 1 /im (see Fig. H) . 

Although there is considerable variation in the quanti- 
tative spectral energy distributions of quasars, the inflexion 
at 1 ^m appears to be a universal feature (Neugebauer et 
al. 1987; Sanders et al. 1989; Elvis et al. 1994), except for 
blazar-type sources where synchrotron emission dominates 
even at optical-infrared wavelengths. The rise to longer 
wavelengths is believed to be due to dust at or near its sub- 
limation temperature (T ~ 1500 K). Kobayashi et al. (1993) 
have successfully modelled the near-infrared continua of low- 
redshift quasars in this manner, and Barvainis (1987) has 
shown that the entire infrared spectrum from 1-100 /im can 
be attributed to dust at a range of temperatures, with the 
hottest dust contributing most strongly at the shortest wave- 
lengths. If the dust is heated by, and in thermal equilibrium 
with, the radiation field of the quasar, then the hottest dust 
must originate close to the central engine. Thus even though 
emission from the hot dust may be optically thin and hence 
isotropic, it could easily be occluded by the larger-scale ob- 
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Figure 2. Comparison of the spectral energy distributions SEDs 
of stellar populations (of various ages) and quasars (subject to 
various amounts of reddening). Stellar population loci (dashed 
curves) are 1-Gyr bursts at ages 1,2,3, and 5 Gyr (younger ages 
correspond to bluer colours) from the GISSEL96 models (Bruzual 
& Chariot 1993,1998). Quasar loci (full lines) have been cal- 
culated using the canonical spectrum described in Section 1 as 
seen through 0,2,7, and 15 mag of visual extinction in the rest- 
frame of the quasar host galaxy (lower extinctions correspond to 
bluer colours). All the SEDs have been normalized at 1 foa. Note 
the similarities between the SEDs of lightly-reddened quasars 
{Ay ~ 2) and old stellar populations blueward of 1 /im, but the 
divergence at longer wavelengths. 



scuring material so that the emergent flux becomes a strong 
function of viewing angle. The attenuation suffered at an- 
gles greater than 6c should, however, be much less severe 
than at shorter wavelengths. The apparent constancy of the 
near-infrared continuum slope in quasars (Neugebauer et al. 
1987), and the possibility that the emitting regions are op- 
tically thin makes it an ideal region of the spectrum from 
which to determine the reddening. 

The structure of this paper is as follows. In Section 2 we 
describe our selection of targets for thermal-infrared imag- 
ing. In Section 3 we describe our new imaging observa- 
tions with the 3.9-m United Kingdom Infrared Telescope 
(UKIRT). In Section 4 we present the results of these obser- 
vations and our analysis of the data. In Section 5 we widen 
our analysis to include archival data from the Hubble Space 
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Telescope (HST), as well as a new near-infrared spectrum 
obtained at UKIRT. In Section 6 we discuss the implications 
of our results, including both a discussion of the effects of 
hidden quasar nuclei on the K-z relation and a test of the 
receding torus model. In Section 7 we review the prospects 
for future observations of the type discussed in this paper. 
Throughout this paper we have assumed Ho = 
50kms~^ Mpc~^, go ~ 0.5 and A = 0. The convention for 
spectral index, a, is that Si, oc v~" , where 5*1^ is the flux den- 
sity at frequency v. We have adopted a canonical rest-frame 
spectrum for a quasar of the following form: the composite 
quasar spectrum from Francis et al. (1991) for wavelengths 
from 0.1-0.6 fj,m, a power law with a = 0.3 from 0.6-1 /im, 
and a power law with a — 1.4 from 1-3 /im. The motivation 
for the broken power-law form at near-IR wavelengths is 
provided by the data of Neugebauer et al. (1987). We adopt 
the empirical extinction curve for the Milky Way given by 
Pei (1992) as our assumed dust extinction curve. 



2 SELECTION OF TARGETS 

Our primary goal in this project was to study a com- 
plete subsample of radio galaxies from the 3CRR sample 
— which we hereafter refer to as the LRL sample (after 
Laing, Riley & Longair 1983) in order to more easily dif- 
ferentiate it from the 3CR sample — in the redshift range 
0.65 < z < 1.20. These objects have 178-MHz radio lumi- 
nosities Li78 > 10^^ WHz~^ sr~^, and they lie at redshifts 
where 3.8 /im imaging probes rest-frame wavelengths in the 
region of primary interest (see Fig. 0). However, a combi- 
nation of scheduling, technical problems, and poor weather 
meant that we were often either unable to observe all the 
sources within a given patch of sky, or conversely were short 
of targets. In the former case, we chose the targets to observe 
based solely on position in the sky so as to make our observ- 
ing programme as efficient as possible, while in the latter 
case, we supplemented our list from the 3CR catalogue of 
Bennett (1962a,b), on which LRL is based. The declination 
limit of UKIRT forced us to exclude targets with 5 > -1-60° . 
Over the course of the three observing runs assigned to 
this project, plus some additional spare time on three other 
nights, we were able to observe all twelve radio galaxies from 
3CR and LRL satisfying a > le*" and a < 04''20™. Our cov- 
erage in the RA range 09*^ < a < M'' is less complete. 
We observed a randomly-selected subsample of five 3C ra- 
dio galaxies from the eight with 0.8 < z < 1.0; the excluded 
sources were 3C 237 (in 3CR only) and 3C 263.1 and 3C 280 
(also in LRL). We provide relevant data for our 3CR (but 
non-LRL) targets in Table hJ; a useful recent reference for 
data on the LRL sample is Blundell et al. (1998). A full list 
of our targets is included in Table H 

In Fig. H, we compare the location of our UKIRT targets 
with the full LRL sample in a plot of radio core luminos- 
ity versus redshift. This figure introduces two new classes 
containing objects traditionally classified as radio galaxies 
(e.g., in LRL) but with broad emission lines in their optical 
and/or near-infrared spectra. The first new class accounts 
for objects with prominent broad lines in their direct-light 
spectra but nuclear optical continuum luminosities below 
the Ms = — 23 limit often taken as the minimum required 
for a quasar (hereafter class-'B' radio galaxies; see Willott et 
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Figure 3. Radio core luminosity (evaluated at or near 5 GHz) 
versus redshift for all objects from the LRL sample with 0.65 < 
z < 1.20 plus additional UKIRT targets. UKIRT targets are in- 
dicated by large symbols: large circles denoting LRL objects; and 
large squares, objects which are only in SCR. The smaller sym- 
bols have meanings as follows: filled squares represent quasars 
(class 'Q'); filled triangles, radio galaxies with broad lines promi- 
nent in direct-light spectra ('B'); filled stars, radio galaxies with 
broad lines seen only because of a large scattered (polarized) com- 
ponent ('S'); and open squares, narrow-line radio galaxies ('G'). 
The numbers in each category in each of three marked redshift 
bins are summarized at the top of the plot, with the first number 
relating to the total number in LRL, the second to the number of 
LRL objects observed with the UKIRT, and the third number to 
the SCR (only) UKIRT targets. A few objects of special interest 
are labelled. 



al. 1998b for further discussion). Since these criteria ignore 
any reddening corrections, class-'B' radio galaxies include 
both lightly-reddened quasars and those that are simply in- 
trinsically weak. The second class accounts for objects with 
broad emission lines which are prominent only in scattered 
light as evidenced by spectropolarimetric observations (here- 
after class-'S' radio galaxies; e.g., Tran et al. 1998). Class- 
's' objects obviously harbour quasar nuclei, but according 
to unified schemes may be being observed very far from a 
pole-on orientation. Fig. M illustrates the dichotomy between 
the radio core properties of quasars and narrow-line radio 
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Table 1. Relevant data on the UKIRT 3CR targets which are not in the LRL sample. Positions are either optical (from McCarthy et 
al. 1997) or radio (from Strom et al. 1990). We provide radio core measurments (and the frequency at which these measurements were 
made) and the largest (projected) angular sizes of the radio sources from the following references: B94, Rogers et al. (1994); L75, Longair 
(1975); R96, Rhee et al. (1996); S90, Strom et al. (1990). 



Name 



Pos 



a (B1950) 5 



Score (mjy) u (GHz) LAS (") Ref 



3C 44 O* 

3C 54 O 

3C 107 O 

3C 114 R 

3C 272 R 



01 28 45.12 +06 08 17.7 

01 52 25.91 +43 31 20.9 

04 09 49.88 -01 07 10.6 

04 17 29.05 +17 46 49.0 

12 22 00.81 +42 23 12.3 



< 1.0 


8.0 


65 


B94 


< 25.0 


5.0 


54 


L75 


< 13.0 


1.4 


13 


R96 


16.0 


5.0 


54 


S90 


0.3 


5.0 


59 


S90 



'Position estimated from the countour plot in McCarthy et al. (1997) since the position given in their Table 3 is incorrect. 



galaxies and hints that class-'B' (and maybe class-'S' 
galaxies are intermediary cases. 



radio 



We next consider possible biases introduced by our in- 
complete UKIRT observational programme. The summary 
numbers included in Fig. H show that we have obtained 
UKIRT data for about half of the possible target LRL ra- 
dio galaxies with no significant redshift bias across the full 
range. There are insufficient data to consider any more than 
one redshift bin in our later analysis, but we note that 
on the basis of Fig. B^ there is very little evidence for any 
change in quasar fraction with redshift across this bin (al- 
though see Lawrence 1991). Again because of lack of data 
we are forced to include the non-LRL objects in our analysis 
without considering them separately from the LRL objects. 
There is some worry that this might introduce a small bias 
because these objects were typically excluded from LRL be- 
cause their accurate radio flux densities fall slightly below 
the LRL limit, and as discussed in Section nl their lower ra- 
dio luminosities might be expected to translate into lower 
quasar luminosities. To assess the impact of this effect we 
plot in Fig. ^the [O II] narrow- line luminosity, L[o n], versus 
178-MIIz radio luminosity, Lits, for the parent LRL sample 
plus our additional UKIRT targets. This plot is relevant 
because according to unified schemes the [O II] luminos- 
ity should be much more closely related to the underlying 
quasar luminosity than Lira (e.g., Hes, Barthel & Fosbury 
1993). Although Simpson (1998) has suggested that the op- 
tical [O III] doublet provides a more reliable probe of the 
intrinsic quasar luminosity, the wavelengths to which these 
lines are redshifted for sources in our sample mean there 
exist very few measurements of their strengths. The effects 
of using different narrow lines are subtle compared to the 
points we shall make on the basis of Fig. W. Any correlation 
in Fig. W is much weaker than the one seen when objects 
spanning several orders of magnitude in Liyg are plotted 
(e.g., Willott et al. 1998c). This demonstrates that the scat- 
ter in the Lyo nj, Lng, relation is sufficiently large that the 
factor ~ 2 difference in Lits between the SCR and LRL 
objects has negligible effect on the expected difference in 
underlying quasar luminosities. Kolmogorov-Smirnov tests 
show that the distribution of I/[o ii] for the UKIRT targets is 
indistinguishable from that for the unobserved LRL galaxies 
alone, and from that for the unobserved LRL galaxies and 
quasars combined. The lower scatter displayed in Lyo n] by 
the quasars is at least in part due to that fact that for a 
significant number of them (7/19) we have inferred the line 
luminosity from broad-band photometry by assuming a sin- 
gle rest-frame equivalent width, whereas in reality there is 



a significant spread in this property. We conclude, and shall 
assume hereafter, that our UKIRT targets form an unbiased 
selection of 3C radio galaxies in the target redshift range. 

In Section owe will need to consider an unbiased sample 
of radio sources including both quasars and radio galaxies. 
To do this we must estimate the number of quasars which 
would exist in a sample with the same incomplete selection 
criteria as our radio galaxy sample. Excluding 3C 454.3, 
whose presence in LRL is due to Doppler-boosted core emis- 
sion leaves 19 LRL quasars with 0.65 < z < 1.20, compared 
to 29 galaxies. This quasar fraction (40 per cent) corresponds 
to a torus opening angle 6^ ~ 53°. Our sample of 19 radio 
galaxies observed with UKIRT (including five from outside 
LRL, but which we have argued are not significantly differ- 
ent) should therefore be accompanied by a 'virtual' sample 
of 12.5 quasars. 



3 IMAGING OBSERVATIONS 

All objects were observed at UKIRT using the IRCAM3 in- 
frared array with a nominal scale of 0'.'281 pixel" ^. Our ob- 
serving strategy consisted of dithering the telescope to place 
the object in nine different locations on the detector ar- 
ray, and after applying a linearity correction to each frame, 
we reduced the data in the following manner. First, each 
of the nine frames was dark-subtracted, and all the frames 
were then scaled to have the same median pixel value and 
median-filtered to produce a flat field. This was normalized 
and then divided into the individual dark-subtracted frames. 
These frames were registered at J and K by determining the 
individual image offsets from the centroids of the brightest 
objects in the fields where possible, and otherwise (e.g. at 
L' , where the high thermal background and smaller field of 
view made it impossible to see anything in the individual 
images) by using the telescope offsets recorded in the im- 
age headers. Each set of nine images was then averaged to- 
gether (ignoring known bad pixels), and objects identified in 
this image as contiguous groups of pixels whose values were 
above a certain threshold determined from the sky level and 
noise. A new flat field was then constructed, ignoring the re- 
gions in the individual frames which contained objects, and 
the above procedure was repeated. Usually we took more 
than one set of nine images, and the final image was pro- 
duced by averaging these sets, after additional registration 
if necessary. Flux calibration was usually performed by ob- 
serving several UKIRT flux standards during the course of 
each night and determining the atmospheric extinction co- 
efficient and array sensitivity in each filter. The L' image of 
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Figure 4. [O II] A3727 narrow-line luminosity versus 178 MHz 
radio luminosity for objects from the LRL sample with 0.65 < 
2 < 1.20 plus the additional UKIRT targets. Symbols have the 
same meaning as in Fig. bl Emission line data are taken from the 
compilation of Jackson & Rawlings (1997) with a few updates 
reported by Willott et al. (1998c); six objects without [O II] mea- 
surements have had their [O III] luminosities scaled according to 
the typical line ratios in McCarthy (1993); an [O II] rest-frame 
equivalent width of 10 A has been assumed for the seven quasars 
without spectrophotometry, and four objects (3C 220.3, 3C 292, 
3C 318 and 3C 325) have not been plotted due to the lack of any 
emission line data. 



3C 22 was calibrated using observations of the single stan- 
dard star HD 225023 taken immediately before and after the 
radio galaxy images. 

Our initial observing run of UT 1995 March 23-26 was 
beset by problems with the IRCAM3 control software, AL- 
ICE. Not only did this result in a significant loss of ob- 
serving time, but our L' images frequently showed residual 
gradients and other structure due to an unstable bias level. 
We excluded from the final coaddition those images which 
were most seriously affected by this problem, but in many 
cases the images could be reduced in the normal manner and 
the only effect was an increased noise level which caused the 



Table 2. Observing log. Listed for each filter and each source are 
the UT date of observation (YYMMDD) and total useful expo- 
sure time (in seconds). /C-band imaging of 3C54 was accidentally 
excluded from the observing programme. 



Name 



J 



K 



L' 



3C 22 


0.935 


970826 


540 


970826 


540 


970130 


540 


3C 34 


0.690 


970825 


1080 


970825 


1080 


970825 


3510 


3C 41 


0.794 


961104 


1080 


961104 


900 


961104 


3240 


3C 44 


0.660 


970824 


1080 


970824 


1080 


970825 


3240 


3C 54 


0.827 


970826 


2160 






970826 


3510 


3C 55 


0.735 


970824 


1080 


970824 


1080 


970824 


3510 


3C 65 


1.176 


961104 


1080 


961104 


1080 


961104 


3240 


3C 107 


0.785 


970824 


1080 


970824 


1080 


970826 


3510 


3C 114 


0.815 


970824 


1080 


970824 


1080 


970825 


3240 


3C 217 


0.898 


950324 


810 


950324 


810 


950324 


1820 


3C 226 


0.818 


950323 


810 


950323 


810 


950323 


1960 


3C 265 


0.811 


950323 


810 


950323 


810 


950324 


2630 


3C 272 


0.944 


950323 


810 


950323 


810 


950324 


2540 


3C 289 


0.967 


950325 


810 


950325 


810 


950325 


2120 


3C 340 


0.775 


970825 


1080 


970825 


1080 


970825 


3510 


3C 352 


0.807 


950325 


810 


950325 


810 


950325 


1980 


3C 356 


1.079 


970825 


4860 


950323 


4860 


970824 


3660 


3C 368 


1.132 


970914 


1080 


970914 


1080 


970914 


3510 


3C 441 


0.707 


970824 


1080 


970824 


1080 


970824 


3390 



overall signal-to-noise ratio to increase more slowly than the 
square root of the total integration time. For this reason, our 
sensitivity at L' during this run was lower than during later 
runs (where this problem was not experienced) , even though 
a number of sources received similar total exposure times. 
We present our complete observing log in Table 0. 

Of those objects observed during the 1995 run, only 
3C 265 was detected, and was extended over several arc- 
seconds east-west. Since this is not a preferred ax;is of the 
host galaxy or radio source, we suspected that the elonga- 
tion may have been caused by telescope drift during the 
two hours over which our observation was made. This was 
confirmed by our 3C 356 images, which were taken over ap- 
proximately two hours of real time, and therefore allowed us 
to measure the true image offsets from sources on the frames 
and compare them to the nominal telescope values. A drift of 
~ 0.9arcsechr~^ in an easterly direction was revealed, and 
the 3C 265 data were re-reduced with this drift accounted 
for. The new image revealed a bright source, this time only 
slightly elongated in an east-west direction. We attribute 
this elongation to imperfect correction of the telescope drift 
and believe the source to be unresolved. We re-reduced the 
other L' datasets incorporating the telescope drift correc- 
tion, but failed to detect any of the remaining objects. The 
addition of the fast guider and tip-tilt secondary mirror for 
the later observations meant that those images were not af- 
fected in this manner. 



4 RESULTS AND ANALYSIS 

In Table we present the results of photometry in 3-arcsec 
and 6-arcsec apertures for the galaxies observed. Although 
some of the limits at L' are fairly bright (L' < 15.5) due 
to the reduced exposure times caused by the problems with 
the instrument electronics, most observations were not af- 
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Table 3. Aperture photometry (converted to zero airmass) from our UKIRT images, 
members of LRL, only SCR. 



All limits are 2.5(t. Sources marked * are not 



Name 



J (3") 



K (3") 



L' (3") 



J (6") 



K (6") 



L' (6") 



3C 22 


0.935 


17.71 ±0.02 


15.74 ±0.03 


13.34 ±0.14 


17.47 ±0.03 


15.57 ±0.03 


13.23 ±0.21 


3C 34 


0.690 


18.53 ±0.02 


16.66 ±0.02 


> 15.72 


18.05 ±0.03 


16.20 ±0.03 


> 14.97 


3C 41 


0.794 


18.38 ±0.06 


16.17 ±0.05 


13.71 ±0.13 


17.93 ±0.06 


15.87 ±0.05 


13.49 ±0.20 


3C 44* 


0.660 


18.24 ±0.03 


16.51 ±0.03 


> 15.75 


17.75 ±0.03 


16.09 ±0.03 


> 15.00 


3C 54* 


0.827 


18.69 ±0.03 




> 15.60 


18.20 ±0.03 




> 14.85 


3C 55 


0.735 


18.92 ±0.04 


17.69 ± 0.04 


> 15.84 


18.57 ±0.05 


17.40 ±0.06 


> 15.09 


3C 65 


1.176 


19.36 ±0.07 


17.25 ±0.05 


15.53 ±0.49 


18.90 ±0.08 


16.83 ±0.06 


> 14.80 


3C 107* 


0.785 


19.24 ±0.05 


17.36 ± 0.04 


> 15.57 


18.95 ±0.07 


17.13 ±0.05 


> 14.82 


3C 114* 


0.815 


19.24 ±0.05 


16.95 ±0.03 


> 15.72 


18.75 ±0.06 


16.46 ±0.03 


> 14.97 


3C 217 


0.898 


19.12 ±0.04 


18.14 ±0.07 


> 15.32 


18.75 ±0.05 


17.68 ±0.09 


> 14.57 


3C 226 


0.818 


18.88 ± 0.05 


17.08 ±0.04 


> 15.41 


18.39 ±0.06 


16.76 ±0.05 


> 14.66 


3C 265 


0.811 


18.09 ± 0.04 


16.67 ±0.03 


14.04 ±0.13 


17.66 ±0.04 


16.35 ±0.04 


13.94 ±0.22 


3C 272* 


0.944 


19.06 ± 0.05 


17.14 ±0.04 


> 15.39 


18.88 ±0.06 


16.86 ±0.05 


> 14.64 


3C 289 


0.967 


18.94 ±0.04 


17.31 ±0.05 


> 15.47 


18.34 ±0.04 


16.79 ±0.06 


> 14.72 


3C 340 


0.775 


18.94 ±0.03 


17.25 ±0.03 


> 15.70 


18.64 ±0.05 


17.12 ±0.05 


> 14.95 


3C 352 


0.807 


18.69 ±0.03 


17.32 ± 0.05 


> 15.41 


18.06 ±0.03 


16.84 ±0.06 


> 14.66 


3C 356a 


1.079 


19.34 ±0.03 


17.57 ±0.03 


> 15.65 


18.92 ±0.03 


17.18 ±0.04 


> 14.90 


3C 356b 


1.079 


19.40 ± 0.03 


17.47 ±0.03 


> 15.65 


18.74 ±0.03 


16.86 ±0.03 


> 14.90 


3C 368t 


1.132 


19.91 ±0.18 


18.12 ±0.15 


> 15.79 


19.26 ±0.29 


17.55 ±0.23 


> 15.04 


3C 441 


0.707 


18.48 ± 0.03 


16.76 ±0.03 


> 15.80 


18.00 ± 0.04 


16.32 ±0.03 


> 15.05 



tWe have subtracted the foreground Galactic star and assumed that the nucleus is at the location identified by Stockton, Ridgway & 

Kellogg (1996). 



fected and the limits provide useful constraints as detailed 
in Section |6[ 

Aperture photometry of our images at J and K agrees 
well with that of Best et al. (1998), with an r.m.s. differ- 
ence of 0.10 mag and no tendency for our magnitudes to be 
systematically fainter or brighter. The agreement with the 
aperture photometry of Lilly & Longair (1984) is also good, 
although the uncertainties in their measurements are often 
large. We note that errors at the level of 0.1 mag are insuffi- 
cient to affect our results. Four radio galaxies were detected 
at L' , and we present our JKL' images of these, together 
with an optical HST image of each, in Fig. y. 

In Fig. y we present our photometry graphically in a 
colour-colour diagram. We also plot the loci for combina- 
tions of reddened quasar plus unreddened starlight. This 
figure can be used to make an approximate determination 
of the nature of individual sources; for example, 3C 22 and 
3C 41 are clearly lightly reddened {Av ~ 2-5) quasars, 
whereas 3C 34 and 3C 44 appear to be normal galaxies 
with no sign of AGN activity in their infrared photometry. 
However, some objects, such as 3C 55, are more ambigu- 
ous, but this ambiguity can be removed by considering both 
the ii'-band morphology and the L'-band magnitude of such 
sources. Considering 3C 55, it is more than two magnitudes 
fainter at K than the locus of 3C quasars (Fig. hi), and it 
is clearly resolved in our near-infrared images, indicating 
that its putative central quasar does not contribute signifi- 
cantly at these wavelengths. The same line of reasoning can 
be applied to the other radio galaxies whose K — L' upper 
hmits do not apparently exclude them from being lightly- 
reddened quasars in Fig. o. The case of 3C 65 is rather more 
ambiguous — for the present, we shall assume that it too is 
lightly- reddened (see also Lacy et al. 1995). 

Although Fig. M can be used to estimate directly the 
extinctions to the quasars in those radio galaxies detected at 



Table 4. Estimated extinctions and unobscured quasar magni- 
tudes, and estimated fraction of non-stellar emission in a 3-arcsec 
aperture at K, determined using the method outlined in Section kl 



Object 



Ko 



L 



K 



3C 22 3 ± 1 

3C 41 6 ± 2 

3C 65 3 ± 1 

3C 265 15 ±4 



14.77 ±0.18 0.90 

15.06 ±0.21 0.51 

17.70 ±0.22 0.28 

14.31 ±0.45 0.14 



+0.10 

0.12 
+0.13 

0.11 
+0.12 

0.09 
+0.05 

0.06 



L' , we opt for a less subjective approach. In this method, we 
adopt various galaxy SEDs (1-Gyr bursts with ages ranging 
from 1 Gyr to the age of the Universe at the source redshift) 
and solve for the quasar and galaxy magnitudes and the 
nuclear extinction, using our JKL' photometry. As Fig. Q 
shows, the J — K and K — L' colours are not especially 
sensitive to age at 2 ~ 1, so there is little variation in the 
derived quantities with the assumed galaxy age. With the 
added constraint that the K magnitude of the host galaxy 
should lie on the K-z relation, we then determine the range 
of each parameter. In Table ]a we present the derived nuclear 
extinctions and unobscured quasar K magnitudes. We also 
provide the estimated fraction of K-h&nd non-stellar light 
in a 3-arcsec aperture. 



5 COMPARISON WITH OTHER DATA 

5.1 Re-analysis of HST data 

Although we have obtained extinction estimates from our 
infrared data alone, they are rather uncertain and we would 
also like to be able to confirm the correctness of our val- 
ues. The easiest way to do this is by extending our study to 
optical wavelengths, and searching for point sources there. 
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Figure 5. Images of the four radio galaxies detected at L' . The grayscale is an optical HST image, and the contour maps are, from 
left to right, J, K, and L', with contours spaced at intervals of 0.25 mag arcsec^^. Each field is 12arcsec on a side. The galaxies and 
the filters in which the HST images were taken are, from top to bottom, 3C 22 (F622W), 3C 41 (F555W), 3C 65 (F675W), and 3C 265 
(F555W). 



However, with the exception of 3C 22, which McCarthy 
(1988) reported to be unresolved in an optical continuum 
image, the quasar nuclei in the other objects are predicted 
to be too faint in the optical {V > 24) to be readily de- 
tectable from ground-based images. We therefore turn to 
the archived Hubble Space Telescope WFPC2 observations 



of these objects (PID 1070, P.I. Longair; see Best, Longair & 
Rottgering 1997) where the point-source sensitivity is much 
greater. The calibrated science data files were retrieved from 
the HST data archive and each pair of images combined and 
cleaned of cosmic rays using the iraf task crrej , with pixels 
unaffected by cosmic rays in both images being assigned an 
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Figure 6. Colour-colour diagram for the galaxies observed, as measured in 3-arcsec apertures, and corrected for Galactic extinction 
using the prescription of Rowan- Robinson et al. (1991). The shaded region on the left indicates the colours of stellar populations formed 
in a 1-Gyr burst with ages 1-5 Gyr, with older populations being redder in J — i^" (as indicated). The solid circle represents the combined 
colours of an unreddened quasar and stellar population. The dashed lines indicate the loci of colours for a reddened quasar and an 
unreddened stellar population of ages 1,2,5 Gyr at z = 1 (there is little variation in the loci with redshift within the range we are 
considering) where the quasar is intrinsically five times brighter than the host galaxy in the X-band. The dotted lines are the loci for 
the case where the quasar is ten times brighter than its host galaxy. Each curve has tick marks at Ay = 2,5,10,15,30 mag. The wedge in 
the top left of the diagram shows the effect of one magnitude of visual extinction at 0.65 < z < 1.20. 



exposure-time-weighted average. Flux calibration was per- 
formed using the on-line WFPC2 exposure time calculator, 
assuming an q = source spectrum. 

We re-analysed these data using our own two- 
dimensional fitting procedure. Heavy smoothing of the data 
revealed that a circularly-symmetric de Vaucouleurs (1948) 
profile provided an adequate fit to the host galaxy compo- 
nent for 3C 22, 3C 41 and 3C 65 while no reasonable galaxy 
model could fit the very peculiar structure of 3C 265 (see 
Fig. H) so no fit was attempted. For the other objects, we fit 
a model WFPC2 point spread function (generated with the 
TINY TIM software; Krist 1995) and a de Vaucouleurs host 
galaxy. The point source flux and the total fiux and effec- 
tive radius of the galaxy component were varied to minimize 
the value of x^- The uncertainty assigned to each pixel was 
calculated as the Poisson error on the source counts added 
in quadrature to the standard deviation of pixel values mea- 
sured in a blank region of the chip. Values of x'^ (and residual 
images) were produced as a function of fitting radius to check 
for systematic effects (such as incorrect background subtrac- 
tion), and as none were found a fitting radius of 15 pixels 



(1.5 arcsec) was used in all cases, and fluxes measured over 
the same aperture. The residual images were structureless 
with the exception of 3C 22 for which the two companion 
objects discussed by Best et al. (1997) were clearly revealed 
(and possibly the F555W image of 3C 41 where there was 
some evidence for excess partially-resolved emission near the 
nucleus) ; in the case of 3C 22 the fitting procedure was re- 
peated after masking out data at the positions of the com- 
panions. The quantitative results of this analysis are pre- 
sented in Table pi also listed are the fluxes estimated using 
the intrinsic nuclear fluxes and extinctions listed in Table W, 
accounting for the Galactic extinctions given in Table H. The 
exclusion of 3C 265 from this analysis is not a great prob- 
lem because its predicted transmitted nuclear flux is several 
orders of magnitude fainter than the detection limit in the 
HST images, and spectropolarimetry (Dey & Spinrad 1996) 
suggests that any compact component is more likely to be 
light scattered close to the nucleus. The agreement between 
predicted and observed HST fluxes is generally very good, 
with the exception of the F555W flux of 3C 41, which is 
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Table 5. The results of our host galaxy plus n ucle ar point source decomposition performed on the archived HST images of the four 
objects detected at L' , as described in Section p.l| . Data obtained on 3C 22 using the F1042M filter have extremely poor sensitivity 
and are ignored. For the best-fitting model in each object /filter combination, we list the reduced x^i the effective radius Tc of the host 
galaxy, the percentage contribution to the HST counts within a 3-arcsec diameter from the galaxy and quasar, and the flux density of 
the quasar component. Note that these percentages typically sum to less than 100 per cent because other components are revealed by 
the modelling process, and evident in the residual map. Significant diffe rences between our best-fit parameters and those derived from 



the same HST dataset by Best et al. (1998) are discussed in Section 
from our near-infrared data (see Section W). 



5.1 



Also listed is the fiux density of the quasar nucleus estimated 



Name 



WFPC2 
Filter 



(arcsec) 



% gal % QSO 



Fluxes (/iJy) 
Observed Estimated 



3C 22 


F622W 


3.3 


0.8 


80 


15 


1.8 


3.0lt»o* 


3C41 


F555W 


2.0 


0.6 


60 


30 


1.3 


O.OS+g:" 




F785LP 


2.5 


0.8 


90 


10 


2.8 


^•09-0.52 


3C 65 


F675W 


1.3 


1.2 


>90 


< 10 


<0.2 


n qQ+O-76 
u.oo_g^23 




F814W 


1.0 


1.5 


>90 


<10 


<0.7 


74+i-2'^ 



predicted to be much fainter than is observed. We discuss 
this diflFerence below. 

It is worth noting some differences in our fitting re- 
sults to those obtained previously by Best et al. (1998) us- 
ing a one-dimensional fitting algorithm on the same dataset. 
While we obtain similar results for 3C 65, our effective radii 
for 3C 22 and 3C 41 are much smaller than those of Best 
et al. We suggest that these authors have overestimated the 
radii in these galaxies because their bright nuclei force them 
to exclude all data within a radius of 1 arcsec, and the fit 
is therefore dominated by the very tenuous detection of the 
low-surface brightness extended starlight — their outermost 
datapoints correspond to fiuxes per pixel well below one- 
tenth of the sky noise. At these levels, one cannot remove 
all the companion objects, and there will also be low-energy 
cosmic rays left in the final image, with the result that the 
surface brightness will be overestimated, leading in turn to 
an overestimate of the effective radius. Fitting the host and 
nucleus simultaneously using the full two-dimensional data 
alleviates all of these concerns, allowing data from the cen- 
tral regions to constrain r^ as well as giving a very low weight 
to individual discrepant pixels at large radii. 

After using the HST data to isolate compact compo- 
nents in the optical it was possible to make more detailed 
studies of the SEDs, and hence better estimates of the nu- 
clear extinctions of our UKIRT L' detections. A minimum- 
X'^ fitting procedure was used to find the most likely values of 
Av, and the fractions of nuclear fiux in 3-arcsec apertures 
at J, H and K. In performing this analysis, we assumed 
that the L' and HST point source fluxes are dominated by 
transmitted light from the reddened quasar nucleus, except 
for 3C 41 where we make a correction for scattered light de- 
scribed below. The results of this analysis are presented in 
Fig. M and Table H. We estimated errors on the fitted quan- 
tities by varying the intrinsic colour of the quasar over the 
range observed in naked quasars accounting properly for the 
photometric error on each (assumed) nuclear flux data point 
(see Fig. m); there were no significant differences in the de- 
rived quantities from fits which used LMC- and SMC-type 
dust extinction curves (from Pei 1992). Note the extremely 
good agreement between the extinction estimates, and espe- 
cially the fractions of non-stellar K-h&nd light, determined 
from the infrared data alone (Table kh and the more refined 
values incorporating the HST data (Table Bl). This indicates 



the ability of infrared imaging alone to provide quantitative, 
as well as qualitative, statements about the nature of radio 
galaxies. Throughout the remainder of this paper, we adopt 
the extinction estimates listed in Table o in preference to 
those of Table [i. Brief notes on an object-by-object basis 
follow. 



5.2 Notes on individual objects 

3C 22 A seemingly straightforward case of a nucleus red- 
dened by Av = 3.6 ± 0.3 which contributes ~ 80 per cent 
of the 3-arcsec aperture K-h&nA light, and dominates the 
unresolved component in the HST data. The lack of optical 
polarization (see Leyshon & Eales 1998) supports this pic- 
ture. The extended near-infrared light is presumably from 
the host galaxy which is clearly visible in the HST image. 

3C 41 As we mentioned in Section ^j|, the observed op- 
tical fluxes of the compact nucleus are much larger than 
we inferred from our infrared imaging, suggesting the pres- 
ence of an additional component, such as scattered light 
from the nucleus. This is supported by the observed po- 
larization properties of 3C 41 reported in Leyshon & Eales 
(1998), where the V-band polarization is P = 9.3±2.3%. For 
electron-scattering over a small range of angles about a mean 
angle, 6, the intrinsic polarization of electron-scattered radi- 
ation is P = 100% X (1-cos^ 6')/(l-|-cos^ 6), and as a lightly- 
reddened quasar, the viewing angle should be only slightly 
larger than the critical angle of the torus, 6c = 53°. There- 
fore the intrinsic polarization of the scattered light should 
be around 50%, and approximately one-flfth of the total 
F555W flux (i.e., ~ 1 fxij) would need to be scattered to 
produce the observed fractional polarization. This fraction 
may be larger if the observed polarization has been diluted 
by virtue of being averaged over a range of scattering angles, 
or if the scattering medium is dust. If the scatterers have an 
asymmetric distribution within the illumination cones, our 
estimate of the scattering angle will be in error, although this 
could serve to either increase or decrease the intrinsic polar- 
ization. In order to perform a quantitative analysis, however, 
we assume that this is the true contribution of scattered 
light. The contribution from scattered light in the F785LP 
fllter is predicted to be similar (from Fig. 4 of Leyshon & 
Eales 1998), and we therefore fit the SED after subtracting 
1 /iJy from each of the two HST point source fluxes to ac- 
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Figure 7. Spectral energy distributions (SEDs-) for 3C 22, 3C 41, and 3C 65. The data have been corrected for extinction in the Milky 
Way using the extinction values given in Table pi and are taken from this paper, Lilly &; Longair (1984) and, in the case of 3C 65, Rigler 
& Lilly (1994). Filled data symbols represent points where we have assumed that all the flux is transmitted light from the quasar nucleus. 
The full line represents a fit using the canonical quasar spectra reddened using a MW-type extinction curve (full lines). The dotted line 
is the most-reddened quasar flt consistent with the data allowing the intririsic SED for the quasar to be bluer by subtracting 0.3 from 
the spectral index of both the power-law components discussed in Section hi, and the dashed line the least-reddened quasar fit consistent 
with the data with a redder intrinsic SED (produced by adding 0.3 to both spectral indices). These assumptions about the spread in the 
intrinsic SEDs are based on the quasar SED study of Neugebauer et al. (1987). 



Table 6. Results of our fits to the observed spectral energy distributions of 3C 22, 3C 41, and 3C 65: the nuclear extinction and fractional 
contributions of the quasar nucleus in the JHK infrared bands in 3-arcsec apertures are listed. The Galactic extinction, Ay (MW), has 
been estimated from the IRAS 100 /.im maps obtained through Sky View (McGlynn, ScoUick & White 1996) and the prescription of 
Rowan-Robinson et al. (1991). Note a typographical error in the value of the Galactic extinction for 30 22 listed in the caption to 
Figure 3 of Rawlings et al. (1995). 



Name Ay (MW) 



fj 



fh 



fK 



30 22 


0.66 


3.6 ±0.3 


0.30 ±0.10 


0.70 ±0.20 


0.85 ±0.15 


30 41 


0.20 


4.4 ±0.4 


0.25 ±0.05 


0.35 ±0.10 


0.75 ± 0.20 


30 65 


0.15 


4.8 ±1.3 


0.06 ±0.05 


o.iot 


0.25 ± 0.05 



IFoUowing Rigler . 



Lilly (1994) we have assumed that 10% of the ff — band flux measured by Lilly . 
component arising from reddened transmitted quasar light. 



Longair (1984) lies in a nuclear 



count for the scattered component. We find Ay ~ 4.4 ± 0.4 
and a transmitted quasar ligiit contribution of ~ 85% to 
the 3-arcsec aperture A'-band emission. As was the case for 
3C 22, the extended near-infrared light can be ascribed to 
the host galaxy mapped in the HST image. 

3C 65 There are two reasons why it remains uncertain 
whether we have detected non-stellar emission at L' in this 
object. First, the detection itself is significant only at the 
2.7a level and then only in the 3-arcsec aperture. Secondly, 
the separation between the location of 3C 65 and possible 
loci of stellar populations in Fig. H is only marginally signifi- 
cant. However, both Lacy et al. (1995) and Stockton, Kellogg 
& Ridgway (1995) have argued previously for the presence of 
a quasar nucleus, and most importantly for our purposes an 
analysis by Rigler & Lilly (1994) of a 0.6-arcsec resolution 
H-hand image points to a nuclear component contributing 
~ 10% of the flux. Our SED fit uses this point, the L' de- 
tection and the HST limits to suggest that there is indeed 
a quasar present which has Ay ~ 4.8 ± 1.3. As suggested 
previously by Lacy et al. (1995), 3C 65 seems seems have an 
unusually faint quasar nucleus for a powerful radio source. 
In the absence of reddening we estimate that it would have 
AIb = —24.8, 1.5 mag fainter than the intrinsic luminosity 
of the 3C 22 quasar and 0.6 mag fainter than the 3C 41 
quasar. Note that 3C 65 is a prominent outlier in Fig. H in 
the sense that its narrow emission line luminosity lies factors 



of 26 and 11 below the line luminosities of 3C 22 and 3C 41, 
respectively. 

3C 265 Although the complex optical structure of 3C 265 
excludes an analysis of the HST data, a visual inspection 
of the images indicates that no point source is present at 
the location of the nucleus, supporting our higher extinc- 
tion estimate of Ay ~ 15 for this source. The fact that 
3C 265 is convincingly detected at L' despite such large ex- 
tinction is undoubtedly due to the very luminous quasar it 
contains, as evidenced by its line luminosity (Fig. M). Al- 
though this source displays strong, extended optical polar- 
ization (Jannuzi & Elston 1991), the near-infrared polariza- 
tion is lower (~ 5% in the H-hand; Jannuzi, private com- 
munication). The thermal-infrared polarizarion is likely to 
be even lower, since the overall SED of 3C 265 rises more 
rapidly in the thermal-IR than a quasar SED. Even as- 
suming electron scattering, the scattered nuclear component 
cannot contribute more than ~ (3//)% of the L'-band flux, 
where / is the fractional polarization of the scattered com- 
ponent, unless it too is heavily reddened. We are therefore 
confldent that our L' detection of this source represents the 
true, transmitted quasar continuum. 
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Table 7. Results of four-Gaussian fit to the continuum- 
subtracted spectrum of 3C 41. The hne widths have not been 
corrected for the instrumental resolution of 300kms~^. 



Line 



Redshift 



Flux 
(10-19 Wm-2) 



FWHM 
(kms-i) 



Han 


0.7929 ± 0.0005 


5.06 ±0.81 


771 ±111 


[N II] -^ 


0.7911 ±0.0008 


2.83 ± 0.74 


695 ± 178 


Hat 


0.7964 ±0.0010 


21.25 ±1.67 


3438 ± 335 



°A6584 only. 

5.3 Ha spectroscopy 

As we mentioned in Section ul, a low value of the nuclear red- 
dening may be indicated by broad wings on the Ha emission 
line, which is redshifted into the J window for 3C 22, 3C 41 
and 3C 265 (but lies where the atmosphere is opaque for 
3C 65). Economou et al. (1995) and Rawlings et al. (1995) 
have both shown the existence of broad Ha in 3C 22, and in- 
ferred low values of the extinction {Ay < 4 and Ay = 2 ± 1, 
respectively), in agreement with our own estimate. Given 
the relatively low line of sight extinction derived for the 
quasar within 3C 41, we took a J-band spectrum with CGS4 
on UKIRT on the night of UT 1997 August 12 to look for 
broad Ha. We used the standard 'ABBA' technique (e.g., 
Eales & Rawlings 1993) for observing the source, keeping 
it on the array at all times so that all 48 minutes of inte- 
gration contributed to the object signal. The 75 lines mm~^ 
grating was used in second order with the short camera and 
a 1-pixel (1.22-arcsec) slit to provide a resolving power of 
A/AA « 1000 at the wavelength of Ha. We used 2x2 over- 
sampling and the spectra were interleaved with bad pixel 
masking using the CGS4DR software package before being 
exported to iraf for coaddition, sky-subtraction, and ex- 
traction of the one-dimensional spectra. The final reduced 
spectrum (extracted from a 1.22-arcsec square aperture) is 
presented in Fig. w. Broad wings are clearly visible on the 
Ha line. The convolution of this spectrum with the UKIRT 
J filter gives a synthetic magnitude J = 19.13, which agrees 
well with the value J = 19.29 determined from our image, 
given that the image was taken in poorer seeing. 

We have analysed the spectrum of 3C 41 to derive an 
independent extinction estimate. First a fit was made to the 
continuum using ten cubic spline pieces and subtracted from 
the data. We then used the SPECFIT package (Kriss 1994) to 
fit the Ha+[N II] blend with a combination of four Gaus- 
sians, representing the three narrow lines and a broad Ha 
component. We fixed the relative wavelengths and intensi- 
ties of the [N II] lines and demanded that they have the same 
width, but allowed the other parameters to vary freely. The 
results are presented in Fig. and Table fl 

The results of our line fitting are very similar to those 
found by Economou et al. (1995) for 3C 22, namely: a rest- 
frame equivalent width of 470 A for the blend (c/. 540 A for 
3C 22); a relative redshift of 1050 kms"^ (960 kms"^) for 
the broad Ha line with respect to the narrow line; and ratios 
[N II] A6584/Ha = 0.56±0.17 (0.59±0.14) and Hab/Ha„ = 
4.19±0.75 (3.92±0.49). Assuming an intrinsic ratio of broad 
to narrow Ha of 40 (e.g., Jackson & Eracleous 1995), this 
implies Av « 3 as the differential extinction between the 
broad and narrow lines, in fair agreement with our more 
refined estimate from the spectral energy distribution. 




Figure 8. J-band spectrum of 3C 41. Note the broad base to the 
Ha+[N II] complex and the apparent detection of [S II]. 




Figure 9. Continuum-subtracted J-band spectrum of 3C 41 and 
best fit of Gaussian lines, as described in the text. 



While Economou et al. (1998) have also used CGS4 to 
detect broad wings on the Ha line in 3C 41, their spectrum 
of 3C 265 does not reveal a broad line. This non-detection 
is in agreement with the much larger extinction we derive 
to the quasar in this object. 



6 DISCUSSION 

In Fig. hd we plot the L' magnitudes and limits of our 
UKIRT targets against redshift, and compare these data 
with reddened quasar spectra. We also plot the estimated L' 
magnitudes of the LRL quasars in the same redshift range, 
extrapolated from optical (usually V^-band) magnitudes, as- 
suming zero reddening. Because the differential reddening 
between (observed-frame) V and L' is large (approximately 
2 mag for Ay ~ 1), even modest amounts of extinction will 
move the quasars up the diagram significantly. Indeed, the 
two apparently faint quasars at high redshift are 3C 212 and 
3C 190, which are both intrinsically red (Smith & Spinrad 
1980) and will therefore be much brighter at L' than we have 
estimated. 

Bearing this in mind, it is clear that our L' limits are 
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Figure 10. The L' versus redshift plane for the LRL quasars plus 
our sample of radio galaxies; symbols have the same meaning as in 
Fig. bI L' magnitudes for the LRL quasars have been estimated 
from the optical magnitudes listed in LRL (except for 3C 190 
where LRL list a photographic magnitude, and we elect to use 
the HST magnitude of de Vries et al. 1997), using the canonical 
quasar spectrum to extrapolate to longer wavelengths. We exclude 
3C 343 from this plot since, although classified as a quasar by 
LRL, it is clearly resolved by de Vries et al. (1997). The solid 
lines show the loci of quasars with the canonical spectrum and no 
reddening for two values of absolute magnitude {Mb = —26 for 
the lower curve, Mg = —23 for the upper). The dashed lines show 
the loci of the same quasars reddened internally by Ay = 2, 7 and 
15. 



much fainter than the loci of unreddened quasars with lumi- 
nosities typical of powerful radio sources at these redshifts. 
We have also been able to detect the weakest quasars asso- 
ciated with such objects, like 3C 65, although at poorer lim- 
iting values of Ay ^ 5. However, the distribution of narrow 
emission line luminosities (Fig. M) suggests that there are 
very few, if any, other objects as weak as 3C 65 in the LRL 
sample at these redshifts. Assuming then that the quasars 
within the many radio galaxies for which we only have limits 
at L' have luminosities typical of those within other power- 
ful radio sources, Fig. [10 makes it clear that the extinction 
must be large in those sources we failed to detect at L' . 

Our data can be used to quantify the size of the dust- 
reddened quasar population. Since all our sources are steep- 
spectrum objects, red objects in our sample are unlikely to 
be due to beamed synchrotron radiation, as has been sug- 
gested for the red flat-spectrum quasar population (Serjeant 
& Rawlings 1996; Benn et al. 1998). Among the 18 sources 
classified as quasars in Fig. [UJ, there are two (3C 190 and 
30 212; Smith & Spinrad 1980) that fall into this category. 
Scaling this to our complete sample of 12.5 'virtual' quasars, 
and adding in 30 22, 30 41, and 30 65, we find that 28li^% 
of quasars are dust-reddened (90% confidence limits) . 



6.1 Influence of quasar nuclei on the K—z relation 
for radio galaxies 

Our first goal for this project was to investigate whether 
the presence of reddened quasar light accounts for a signifi- 
cant fraction of the total K-h&nd flux in 30 radio galaxies. 
This could provide an explanation for the brighter K mag- 
nitudes displayed by 30 radio galaxies compared to those 
from the fainter 60 and B2 samples (Eales et al. 1997). We 
have therefore undertaken a similar analysis to that of Eales 
et al., which involves fitting the least-squares regression line 
to the 30 AT-band data, and then examining the residuals 
about this line for both the 30 and 60 samples (we do not 
supplement the fainter sample with galaxies from the B2 
survey). We use the same 63.9 kpc metric aperture as Eales 
et al. and use the same prescription for correcting measure- 
ments made in apertures of other sizes. 

We first confirm the representative nature of our sample 
by fitting a regression line to the if-band magnitudes of our 
18 galaxies. All 13 galax;ies from the 60 survey in the same 
redshift range (0.65 < z < 1.20) lie above this line (i.e., they 
are fainter), with the mean difference being 0.72 mag. This is 
in good agreement with the 0.59 mag discrepancy found by 
Eales et al. (1997), and we can use a Kolmogorov-Smirnov 
test to rule out the hypothesis that both samples are drawn 
from the same distribution at better than 95% confidence. 

We then correct the photometry for the presence of red- 
dened nuclear light, using the fractional contamination de- 
termined in Section 5.1 and repeat the process. This time 



^^ Wm ^ to provide 
this is several times 



the 60 galaxies are on average only 0.56 mag (40 per cent) 
fainter, although the samples are still different at greater 
than 90% confidence. Our confidence levels are lower than 
those of Eales et al. (1997) solely because our samples each 
contain only half the number of objects as theirs. 

Ooupled with the conclusion of Leyshon & Eales (1998) 
that scattered light is not a major contributor to the to- 
tal A'-band flux of most 30 radio galaxies, our analysis is 
able to rule out quasar light as the cause of the correlation 
between radio luminosity and near-infrared magnitude at 
z ~ 1. Emission lines can also be ruled out since they would 
need to have a total flux of ~ 2 x 10 
one-third of the total -R'-band light — 
greater than the flux of the [O III] AA4959,5007 doublet yet 
the strongest lines which fall into the A'-band over our red- 
shift range are intrinsically much fainter (see also Rawlings 
et al. 1998a). Finally, the nebular emission which has been 
suggested provides a significant fraction of the optical light 
and might explain the alignment effect (Dickson et al. 1995) 
is a negligible contributor at near-infrared wavelengths. As- 
suming [O Il]/H/3 -^ 4 (McOarthy 1993), it can contribute 
no more than about one per cent of the total A'-band flux. 

The 'extra' component in 30 radio galaxies, compared 
to 60 radio galaxies, is therefore not directly related to the 
more luminous central quasar, and must be additional stellar 
luminosity, and presumably stellar mass, in the host galaxy. 
This fits together with the smaller linear sizes of 60 radio 
galaxy hosts measured by Roche et al. (1998). Plausible rea- 
sons why a correlation between host galaxy mass and radio 
luminosity is seen ai z ^ 1, but not at lower redshifts, are 
given by Best et al. (1998). 

Finally, we note that the r.m.s. deviation of the 30 
galaxies about the regression line is reduced from 0.59 to 
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Figure 11. Modified A'-band Hubble diagram for radio galax- 
ies from the LRL and 6C samples (open circles and crosses re- 
spectively). The solid circles are our sample of 18 galaxies with 
X-band photometry, after correction for reddened quasar light. 
The two vertical lines delineate the redshift range of our sam- 
ple (0.65 < z < 1.20) and the heavy line between them is the 
least-squares fit to the corrected magnitudes. See Figure hi for 
references for the photometry. 



0.43 magnitudes when the correction for non-stellar contam- 
ination is made. This reduced scatter is the same as that ob- 
served around Bales et al.'s (1997) least-squares regression 
line for z < 0.6. 



6.2 Test of the receding torus model 

In order to be able to compare the receding torus model 
to our observational data, we have to associate lower lim- 
its to the nuclear extinction with the radio galaxies which 
we failed to detect at L' . Estimating these extinctions is a 
very uncertain affair, since we do not know their intrinsic 
luminosities. Table indicates that a 'typical' source might 
have L'q ~ 13 {K — L' — 1.87 for our canonical spectrum 
at all redshifts z < 1.2), yet the quasar within 3C 65 ap- 
pears to be nearly ten times fainter, and the z — 0.768 
quasar 3C 175 is a magnitude brighter (Simpson, unpub- 
lished). We adopt Lq = 13.5 so that the undetected sources 
must have Ai^i > 2, or Av ^ 15 (at z = 1.2; at lower red- 



shifts, the rest-frame visual extinction will be higher for the 
same observed-frame yl^,/) and use this limit for all 15 of 
our non-detections. Fig. nn clearly shows the conservative 
nature of this assumption. 

The mean radio luminosity of the radio sources with 
0.65 < z < 1.20 plotted in Fig. W is about 70 times greater 
than that of Hill et al.'s (1996) 0.1 < 2 < 0.2 sample. HiU 
et al. infer an opening angle of 20° for their sample and, 
applying the relationship tan^c oc Liyg, we expect dc = 52°, 
in very good agreement with the value derived directly from 
the quasar fraction in LRL. Using the same value for the 
density of the torus, nn, as Hill et al., the receding torus 
model predicts the distribution of extinctions shown by the 
dotted line in Fig. nd. Although the model does predict fewer 
lightly-reddened sources than at low redshift (~ 30% with 
Av < 6 mag, as opposed to ~ 50%), it still overpredicts 
the number. By increasing the density of the torus by about 
an order of magnitude above that derived for Hill et al.'s 
sample, we can fit the data of Figure ha with the dashed line. 
Since this model is required to reproduce both the numbers 
of naked and lightly-reddened quasars, there is very little 
additional information in Fig. h2^ with which to compare 
its predictions. We cannot therefore claim any quantitative 
agreement between the model and our observations. 

There are several points to be made about the receding 
torus model. First, one nust be careful when considering 
objects with low extinction, since there may be additional 
obscuration which should not be attributed to the torus. The 
line of sight to the quasar nucleus could easily suffer a few 
magnitudes of visual extinction in its passage through the 
host galaxy, or through foreground companion objects such 
as those seen in our images of 3C 22. The one-sided radio 
jet displayed by 3C 22 suggests a smaller angle between the 
radio (and torus) axis and the line of sight than in 3C 41 and 
3C 65, which have similar nuclear extinctions, so it is quite 
possible that the torus provides only a fraction (if any) of 
the obscuration. In fact, of the four quasar nuclei we detect 
at L' , it is only 3C 265 that we can say with reasonable 
certainty is seen through the heavy obscuration of the torus. 

The fact that we have needed to increase one of the 
parameters of the receding torus model, namely the torus 
density nn, to allow the model to fit our data, suggests that 
some refinements are needed. The torus density could in- 
crease as a function of redshift, or quasar luminosity (either 
directly or through the mass of the central black hole or 
host galaxy), or both. Obviously this question can only be 
answered by examining the distribution of extinctions in a 
sample selected at a fainter radio fiux limit, such as the 
7C sample (selected at a flux level 25 times fainter than 
LRL; Willott et al. 1998a) although such a study is clearly 
beyond the capabilities of UKIRT. We can, however, glean 
something just from the quasar fraction in the 7C sample in 
the same redshift range, which is consistent with the model's 
prediction of 11% (Willott et al. 1998a). The relatively smaU 
number statistics prevent this from being a definitive test, 
although the agreement is certainly encouraging and sug- 
gests that the torus height might be the same in all sources, 
irrespective of their redshifts and luminosities. 
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Figure 12. Cumulative distribution of nuclear extinction for our 
sample of radio sources (including the 12.5 'virtual' quasars). The 
solid line shows the data and the dotted line shows the prediction 
of the receding torus model of Hill et al. (1996). The dashed line 
shows the prediction for the receding torus model, renormalized 
to the statistics of our sample. 



7 UTILITY OF THE THERMAL IMAGING 
TECHNIQUE 

The sensitivity of thermal imaging observations to point 
sources is essentially limited by two factors. First, the size 
of the detection aperture controls the detection threshold, 
and so the better the image quality achievable, the greater 
the sensitivity. Secondly, there is effectively a maximum in- 
tegration time set by the stability of the background from 
the telescope, and so the rate at which source photons are 
detected is also important. Obviously then, the diffraction- 
limited images that the new class of 8-m telescopes should 
be able to obtain at L' (FWHM < 0'.'2) will enable much 
deeper observations than we have been able to undertake. 

The Infrared Camera and Spectrograph (IRCS; Toku- 
naga et al. 1998) to be installed on the 8.3-m Subaru Tele- 
scope in 1999 will be able to achieve a 3(t detection of a 
point source with L' = 19.0 in 1 hour of integration, by us- 
ing a 0'.'5 detection aperture. For an intrinsic magnitude of 
Lf) « 13, this would allow the detection of the quasar nu- 
cleus through as much as Ay ~ 50 mag of obscuration. In 
reality, the sensitivity will be slightly worse than this due to 
the presence of the host galaxy which is expected to have 



L' ~ 17-18 within a O'.'S aperture (assuming r^ = I'.'S and 
K — L' = 1). However, the excellent image quality should 
enable the nuclear point source, if present, to be separated 
from the underlying galaxy with relative ease and so much 
of the increased sensitivity will be retained. 

The ability to measure obscuring columns to this depth 
would provide a much more rigorous test of the receding 
torus model than could be presented here. The model shown 
by the dashed line in Fig. [12 predicts only six sources should 
be detected with 15 < Ay < 50, while the remaining nine 
will hav e even larger obscurations. As was discussed in Sec- 
tion 6.2, it is important not to put too much weight on those 



sources with low extinction {Ay ^ 5) since some or even all 
of this material may come from a source other than the 
torus. Given the apparent evidence that no property of the 
torus is constant, the predictive power of the receding torus 
model rests in its ability to infer the obscuration towards 
the nuclei which are too heavily obscured to be detected 
directly. 

Improved sensitivity of the thermal-infrared imaging 
method will also allow studies at much higher redshifts. This 
is an important extension because there are obvious difficul- 
ties with forming enough dust for an Ay ~ 50 torus in an 
arbitrarily short time, so there is likely to be some critical 
redshift beyond which all radio sources have low Ay. 

Thermal-infrared imaging may also play an important 
role in understanding the intrinsic orientation dependence 
of quasar emissions. If the thermal-infrared flux arises from 
optically-thin material and is isotropic, then statistical stud- 
ies of, for example, the B — L' colours of z ~ 1 quasar nuclei 
could be used to test whether the (rest-frame) UV contin- 
uum is also isotropic. Low-frequency selected samples, such 
as those from the 3C and 7C samples, should contain quasars 
randomly distributed in orientation with 9 < Oc, and would 
provide a suitable basis for such a study. One possibility 
would be to look for a correlation between B — L' colour 
(corrected for reddening) and radio core dominance. Vari- 
ous studies (e.g., Netzer 1987) have suggested that the UV 
continua, and the low-ionisation broad emission lines arise 
in optically-thick disks and thus emit anisotropically. 



8 SUMMARY 

We have presented the results of a near- and thermal- 
infrared imaging survey of a representative sample of 19 
z ~ 1 radio galaxies from the 3CR and 3CRR (LRL) sam- 
ples. Four objects were detected at L' , with red K ~ L' 
colours which we attributed to partially-obscured quasar 
light. We have used our infrared data to infer the nuclear 
extinctions and contributions from non-stellar light at K 
for these four sources, and find the results to be in excel- 
lent agreement with a more detailed analysis incorporating 
HST data. We deduce that the fraction of dust-reddened 
quasars at z ~ 1 is 28tll% (90% confidence). We find that 
the non-stellar light provides as much as 80 per cent of the 
small-aperture K-hand flux in two objects classified as 'radio 
galaxies', although it is unable to account for the brighter 
iC-niagnitudes of 3C radio galaxies compared to 6C radio 
galaxies at the same redshift, which must therefore be due to 
a correlation between radio and host galaxy luminosities. We 
have also investigated the observed distribution of nuclear 
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extinctions in our sample in the context of the receding torus 
model. Although the smaller fraction of lightly-reddened nu- 
clei is in qualitative agreement with the predictions of the 
model, a quantitative match can only by made by altering 
one of its parameters. We therefore conclude that the reced- 
ing torus model has limited predictive power when compar- 
ing samples of radio sources at different redshifts and/or ra- 
dio luminosities, although we have been unable to determine 
whether it can be applied to a single sample of sources with 
similar properties. Deeper thermal-infrared imaging, which 
will be possible in the near future, will be able to answer 
this and other questions. 
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